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The study and control of coherence in quantum systems is one of the most exciting recent develop- 
ments in physics. Quantum coherence plays a crucial role in emerging quantum technologies as well 
as fundamental experiments. A major obstacle to the utilization of quantum effects is decoherence, 
primarily in the form of dephasing that destroys quantum coherence, and leads to effective classical 
behaviour. We show that there are universal relationships governing dephasing, which constrain the 
relative rates at which quantum correlations can disappear. These effectively lead to speed limits 
which become especially important in multi-partite systems. 



One of the principle distinguishing features between 
classical systems and quantum systems is the existence 
of quantum coherences leading to correlations that can- 
not be accounted for classically. For example, the phe- 
nomenon of entanglement [JJ and the violation of non- 
local realism [2] are such consequences. The manipula- 
tion and preservation of such coherences is vital for tasks 
such as the construction of quantum information process- 
ing devices [3], quantum communication [5], cryptogra- 
phy [5], and metrology [BJ. It has also been suggested 
that quantum coherence plays a role in certain biological 
processes [7]. Unfortunately, the inevitable interaction 
of quantum systems with the environment leads to deco- 
herence, the dominant form of which is dephasing, or the 
disappearance of the off-diagonal elements of the density 
operator of the system [8]. The rates at which these co- 
herences decay are crucial as they determine how quickly 
the system approaches classicality, a state where quan- 
tum correlations are lost. Considerable effort has been 
expended on understanding the fundamentals of decoher- 
ence and quantum correlations, and how to protect the 
latter and prevent rapid decay using encodings in pro- 
tected subspaces or subsystems, for example [9 . 

Previously, some surprising relationships between the 
dephasing rates in multi-level quantum systems have 
been uncovered [T01 ITT) . These stem from the need to 
preserve positivity of the density operator, or more gen- 
erally the complete positivity of the evolution [JJ], and 
lead to constraints on the relative rates of dephasing. The 
surface of this phenomenon and has only been scratched, 
and in this work we elucidate the universal nature of 
these constraints and present a general framework. An 
important consequence of these constraints are effective 
speed limits on the decay of correlations and entangle- 
ment in multi-partite systems. These speed limits are 
independent of the details of the Hamiltonian evolution 



'Electronic address: daniel.oi@strath.ac.uk 

^Electronic address: sgs29SScam.ac.uk, sgs29@swan.ac.uk 



or decoherence mechanisms and therefore apply to a large 
class of quantum systems from nuclear spins, to atoms, 
ions and quantum dots to biomolecules. 

We consider an iV-level quantum system subject to 
Markovian dephasing whose evolution can be described 
by a Lindblad master equation [131 114] . A pure dephas- 
ing process leaves the populations of the basis states 
invariant, and leads to decay of the magnitude of the 
off-diagonal elements (coherences), as well as frequency 
shifts. Previous work for three-level systems found 
that the decay rates and frequency shifts were con- 
strained [TOJITT]. Additionally, (TU] also gave some partial 
results for four-level systems and showed such constraints 
must exist for higher-level systems but the equations were 
intractable in general. 

Here, we present a canonical form for the pure dephas- 
ing Lindblad operators which allows us to derive a gen- 
eral form for the constraints for iV-level systems. These 
form a hierachy of inequalities, defining a convex cone of 
allowed dephasing rates. The general form also allows us 
to invert physically observed dephasing rates to define a 
unique set of canonical dephasing operators, which reflect 
correlations between noise processes such as fluctuations 
in the energy levels, and may serve as a useful diagnostic 
tool. In multi-partite systems these constraints induce 
speed limits on the decay of non-local quantum correla- 
tions and entanglement in terms of the local dephasing 
rates. 



Results 

Canonical Dephasing Operators 

The key idea is that pure dephasing of an iV-level 
system may be modelled by a diagonal Hamiltonian 
H = diag(A„) and a canoncial set of N — 1 or fewer di- 
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agonal Lindblad operators {Vfe}t =1 of the special form 



Vk = diag ... 0, a 



(fe) 
fe+i> 



,(fc) 

A' 



where the non-zero diagonal elements a n can be com- 
plex except for the hrst non-zero element, a^} v which 
is set to be real and non-negative. The density operator 
elements evolve as 
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with effective frequencies given by u> mn — A m — A„ + 
Aoj mn , and dephasing- induced frequency shifts and de- 
coherence rates 
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Im(a mk a* nk ), 
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(|a m fe| 2 + \a n k\ 2 ) ~ Re(a mfc a* fe ). (3b) 



The populations are constant as uir, 



0. The 



off-diagonal elements decay with the damping rate T mn 

\Pmn(t)\ =e- tr ™"|p mn (0)|. (4) 



If the 

flnfc are real then the expressions simplify, r„ 



and there are no frequency shifts, 
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Auj mn = 0. 

As shown in the Methods, any set of pure dephasing 
Lindblad operators can be transformed to this form leav- 
ing the total superoperator unchanged. This reduces an 
arbitrary number of parameters, specified by the non- 
zero elements of an arbitrary set of dephasing operators, 
to N(N — l)/2 parameters in the canonical form. The 
number of free parameters matches exactly the number 
of dephasing rates T mn and frequency shifts Auj mn for 
an iV-level system. 



Inverting Dephasing Rates 

Using this canonical form, we can determine a set of 
standard operators which generate experimentally ob- 
served dephasing rates T mn and frequencies u; m „(or fre- 
quency shifts). The inversion process relies on the fact 
that the dephasing rates involving the first k + 1 levels 
depend only on the first k dephasing operators, i.e., T 12 
determines the first non-zero element of Vi, which to- 
gether with {T13, T23, AUJ23} then determines a further 
three real parameters, and so forth. Hence, it is a sim- 
ple matter of iteratively solving a nested set of quadratic 
equations, as detailed in the Methods. 

Naturally arising in the course of this inversion, a set of 
constraints on the allowed dephasing rates and frequency 
shifts takes the form of N — 1 inequalities involving the 
first N levels 
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FIG. 1: Convex cone of allowed dephasing rates for N = 
3 and real dephasing, i.e., Ao>23 = 0. The axes are 

x = (r 12 + r 23 + ri 3 )/\/3, y = ^(iha - r 23 ) and % = 

— (ri3 + 1^23) \ and the constraint equation be- 
comes x 2 /2 > y 2 + z 2 , which defines a circular cone in the 
positive octant of in the parameter space of {Fi2, F23, Ti3}, 
tangential to the F12T23, I^Fis, and Fi3Fi2 planes. 

(n) 

where the a\ can be expanded in terms of the T mn i with 
m,n' < n. A more symmetric form of the constraints is 



possible, e.g., for N = 3 

2(ri2T23 + r 2 3ri3 + r^Fis 
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which is reduces to Eq. (25) in [10] if AW23 = 0. However, 
as the number of decoherence rates grows as N(N — 1) /2 
and the inequalities involve products of (N — 1) r mn , 
there is a combinatorial explosion in the number of terms 
in the constraints, which is why previous attempts to ob- 
tain a general form for the constraints failed. For exam- 
ple, the four-level constraint contains 22 terms and the 
five-level constraint contains 130 terms. 

These inequalities form a convex cone of allowed de- 
phasing rates whose boundary is formed by "hypersur- 

faces" defined by — for some k > 1. For N = 3 
there is only a single constraint equation and the con- 
vex cone of allowed dephasing rates can be visualized as 
shown in Figure [T] 



Speed Limits for Entanglement Decay 

The constraints for the decoherence rates and fre- 
quency shifts have important implications for a wide 
range of physical, chemical and biological systems where 
phase relaxation is a dominant process. A consequence 
is the imposition of relative speed limits on the rates at 
which coherences can decay, especially in multi-partite 
systems where entanglement decay is strictly bounded 
above by the single qubit dephasing rates. Dephasing 
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FIG. 2: Constraint violation leads to non-physical states. A 
plot of the minimum eigenvalue of p(t) starting with po = 
!*o}(*o| subject to pure dephasing (H = 0) with F23 = Ti4 = 
pT for different values of p shows the emergence of negative 
eigenvalues for p > 2. 



can be spatially correlated, the Markovian condition only 
constrains the temporal correlations in the noise. In gen- 
eral, the dephasing rates can be of a non-local form. 

Let us start with two qubits where we label the basis 
states by |1) = |00), |2) = |01), |3> = |10) and |4> = |11). 
Assuming that both qubits have the same local dephasing 
rate, i.e. T — T12 = T13 = r 2 4 = 1^34, then the allowed 
decoherence rates for the non-local coherences Ti4 and 
r 2 3 are determined by T. The first non-trivial constraint 
(a 2 3 ^) 2 > gives < r 23 < 4r. The second constraint 

(4 4) ) 2 > leads t0 r 23 + r i4 < 4r [24 . Thus, to en- 
sure complete positivity of the evolution, the non-local 
coherences P23 and P14 can decay at most four times as 
fast as the local coherences, and the sum of the non-local 
decay rates can be no more than 4I\ If they are equal 
r 23 = T 14 = T e we obtain T e < 2T, and Fig. ^demon- 
strates that violation of the bound leads to violations of 
positivity, i.e., non-physical states. 

The limits to non-local coherence decay relate to the 
entanglement between qubits. Starting with the maxi- 
mally entangled Bell state |* ) = ^(l 00 ) + I 11 )), the 
state evolving under pure dephasing 
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has concurrence C(t) = e~ tL 14 [TS], thus T u + T 23 < 4r 
implies that the concurrence cannot decay faster than 
four times the local decoherence rate T. Here, the decay 
of non-local coherences is not lower bounded, i.e. non- 
local coherences can survive indefinitely even for finite 
local decay rates and in this case the entanglement is 



preserved if Ti4 = 0, i.e. there is no sudden death of 
entanglement |17) . 

Alternatively, starting with the maximally entangled 
two-qubit cluster state |<F CS ) = |(|00)+|01)+|10> — 111>), 
we obtain 
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In this case the entanglement can decay even if both 
non-local dephasing rates vanish, Ti4 = T23 = 0, in 
which case the concurrence satisfies 2Ci(i) = |e~ rt + 
1| — |e~ r * — 1|, which tends to zero as t — > 00. If one 
of the two non-local concurrences is 4r and the other 
is 0, e.g., Ti4 = 4r, T 2 3 = 0, the concurrence simi- 
larly decays asympotically but faster. When the non- 
local coherences decay at the same rate Ti4 = T 2 3 = 2r 
we have C 2 (i) = max{0, \{2e- Tt + e~ 2rt - 1)}, and the 
concurrence vanishes when 2e~ rt + e~ 2rt — 1 = 0, i.e. 

= -r^ 1 log(V^ - 1) « 0.383F- 1 , i.e., we observe sud- 
den death of entanglement rather than asymptotic decay. 

We can extend this to n-qubit systems. Given n qubits 
with the same local dephasing rate T, we can simply ap- 
ply the results above to any subsystem consisting of two 
qubits, i.e., the entanglement between any two qubits in 
the system cannot decay faster than 2r. For larger sys- 
tems there are more constraints so in practice the rate 
of entanglement decay between any two qubits would be 
even more restricted. For example for a three-qubit sys- 
tem with have = 2 3 = 8. Assuming the local de- 
phasing rate for each qubit is T, and the dephasing rate 
involving two- and three-qubit transition terms are p\T 
and /x 2 r, respectively, there are 8 constraints restricting 
the allowed values for p,\ and /i 2 . From the constraints 
for the two-qubit system we know that < p\ < 2, but 
Fig. [3] shows that the set of (/ii,/i 2 ) that satisfy all the 
constraints is much smaller. Each additional constraint 
reduces the set of allowed dephasing rates. 



Discussion 

The underlying basis for the dephasing constraints is 
correlations between noise acting on different energy lev- 
els of the system. The canonical dephasing operators 
reflect underlying physical processes with different corre- 
lation properties. For example, a canonical dephasing op- 
erator with a single non-zero element can be interpreted 
as the result of the fluctuation of a single energy level. 
Multiple non-zero diagonal entries correspond to corre- 
lated perturbation of more than one level. These fluctu- 
ations can be correlated across levels even for Markovian 
dynamics where the noise statistics are temporally un- 
corrected. 

An example of where noise correlation can occur is 
magnetic field fluctuations acting on a spin-1 particle 
where the coupling is of the form BzSz- This leads 
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creases, the decay of the non-local coherences becomes 
constrained even more by the local dephasing rates. This 
reflects the general robustness of the non-local correla- 
tions in multi-partite systems [5T] . Conversely, there are 
suggestions that dephasing can play a positive role in bi- 
ological processes [22l E2] • Dephasing has been mooted 
to enhance the transport of energy networks such as pho- 
tosynthetic harvesting complexes. In such systems, mea- 
surement and analysis of the dephasing may illuminate 
these processes and lead to better energy collection de- 
vices. 



O.Z 0.4 0.6 0.8 1 1,2 1.4 1.6 1,8 2 

FIG. 3: Constraint violation map for three qubits with local 
dephasing rates F, two-qubit dephasing rates fJ.iT, and three- 
qubit dephasing rates /^T. For < < 2 the first four of 
the eight constraints are satisfied but additional constraints 
may be violated, e.g. in the yellow region, constraint 5 is 
violated. Each subsequent constraint further restricts the set 
of allowed rates. 



to anti-phase perturbations of the Sz = ±1 levels and 
the canonical operator in the basis {|0),|1),| — 1)} is 
oc diag(0, 1, — 1). If the coupling was instead of the 
form BzS z , then the canonical operator would be oc 
diag(0,l,l). 

Dephasing not only can lead to exponential damping 
of the coherences but also can produce shifts in their fre- 
quencies. Not all of these frequency perturbations can be 
accommodated by modifying the system Hamiltonian in 
general, the residual shifts are intrisic to the decoherence 
processes. Whilst pure damping can be generated by 
phase diffusion due to random drift of the energy levels 
(Wiener-Levy process) |I8j . the frequency shifts can be 
caused by phase kicks, or discrete random phase jumps 
with a poissonian arrival time, though this also produces 
additional damping 19J. Phase kicks occur, for example, 
by collisional processes in gases, whereas phase diffusion 
can be generated by white noise acting on the energy 
levels. 

It is possible to derive dephasing constraints from 
physical models of the noise directly, though deriving the 
general multi-level constraints is considerably more dif- 
ficult than the methods shown here. However, once the 
observed dephasing rates have beeen decomposed into 
their corresponding canonical set of dephasing operators, 
we can assign physical mechanisms by which they occur, 
and hence perform system diagnostis or analysis. The 
ability to identify sources of dephasing will be vital in 
producing coherent quantum devices and improving their 
performance. 

In the context of multi-partite systems, the constraints 
we have derived have implications for the preservation 
of non-local correlations. As the number of parties in- 



Methods 

Canonical dephasing operators. We start with the 
Lindblad master equation (LME) for Markovian open 
quantum system evolution 

p(t) = -i{Hp(t) - p{t)H) + £ D ( P (t)), (6) 

where p(t) is the density operator describing the system 
state (defined on the system Hilbert space H s ), H is an 
effective Hamiltonian H and the superoperator £d{p) 
takes the form £ D (p) = £ fc V[V k ]p with [TJ [H] 

V[V k ] P = V k pVl - \(Vlv kP + pVlv k ). (7) 

for operators 14 onHs- 

The operators (H, V k ) define a pure dephasing process 
with respect to the basis 03 = {1".)}^! if and only if H 
and all V k are simultaneously diagonal with respect to 
03, i.e., we have 

H = J2 A « n « = diag(A„), A„ G R (8a) 

n 

Vfc = ^2 7rifcII rl = diag(7„ fc ). (8b) 

n 

This is easy to see since by definition of a pure dephasing 
process the populations of the basis states remain con- 
stant, and thus each basis state \n) is a steady state of the 
system. This is possible only if the subspace spanned by 
each basis state |n) is Vfc-invariant for all V k [2(3 . This 
shows that all V k must be diagonal in the chosen ba- 
sis. Since II„ is diagonal and diagonal operators com- 
mute we have 2?[Vfc](II n ) = for all n and all k. As 
\n) is a steady state, i.e., IT T1 (t) = 0, it also follows that 
— iHH n + iH n H = for all n. Inserting this into the 
general form of the LME ^ gives the explicit equation 

Pmn(^) y^mn ~t~ fmn) Pmn(t): (9) 

for the evolution of the matrix elements p mn = (m\p\n) 
of the density operator, or in integral form ^ with fre- 
quencies U) mn = A m — A„ + Aw mn and dephasing induced 
frequency shifts and decoherence rates given by ([3|. 

Any set of diagonal Lindblad operators {V k } generates 
pure dephasing dynamics but the set of Lindblad oper- 
ators {Vfc} generating a certain dynamical evolution is 



5 



not unique. In particular, we have unitary invariance, 
i.e., given any set of Lindblad operators {14}, the set of 
operators {Wj} defined by 



(10) 



where Uj k are elements of a unitary matrix, generates the 
same dynamics as X^^I^fc]/ 9 = Sfc^I^fe]/ 9 - Further- 
more adding multiples of the identity, V k — > Vk + oil, to 
a Lindblad operator Vf. only changes the effective Hamil- 
tonian 

V[V k + al]p = \[aVl - a*V k , P {t)]+V[V k ]p, (11) 

and thus the dynamics is unchanged if we replace V k by 
V k + al and H by H + \{aVl - a*V k ). 

The invariance of the dephasing dynamics under the 
two "gauge transformations" (10) and (11) allows us to 



transform any set of dephasing operators {Vi} into an 
equivalent set of dephasing operators in canonical form 
defined in Eq. ([!]), which yield the same observable de- 
phasing rates {r m „} and dephasing shifts {Auj mn }, using 
Algorithm [T] The process is constructive and, using a nk 
intead of a„ , the key steps can be described as follows: 
(1) Using ( [TT| ) we ensure that a\ k = for all V k , mod- 
ifying the Hamiltonian by 



(12) 



as necessary. 

(2) We replace the Lindblad operators V\ = 
diag(0, a 2 i,a 31 , . . .) and V 2 = diag(0, a 22 , b 32 , . . .) with 
a 21 = r 2 ie^ 21 and a 22 = r 22 e^ 22 by {I-Fi, VF 2 } with 

W x = UiiVi + u 12 V 2 = diag(0, c, . . .) (13a) 
W 2 = u 21 V x + u 22 V 2 = diag(0, 0, *, . . .) (13b) 

with the unitary coefficient matrix 



u = 



1 ( r 2 ie l 



c \r 22 e~ 



r 22 e 
-r 21 e 



,#21 
-#22 



(14) 



and c = yj ri> A + r| 2 , which is dynamically equivalent to 
{Vi,V 2 } due to ([Tof. 

This result allows us to reduce an arbitrary number 
of parameters, specified by the non-zero elements of a 
general set {V k } of dephasing operators to N(N — l)/2 
parameters in the canonical form. Note that the num- 
ber of free parameters matches the number of dephas- 
ing rates L mn for an TV-level system. The procedure 
will work, i.e., to produce a set of canonical dephasing 
operators that reproduce the observed dephasing rates 
and shifts, provided that these satisfy the positivity con- 
straints. Furthermore, if the observed dephasing rates 



and shifts lead to constraint violations these will be de- 
tected and flagged, and this information can be used to 
further investigate if the violations can be explained in 
terms of uncertainty in the observed data, e.g., due to 
measurement errors, or if they are indicative of processes 
that would invalidate the Markovian dephasing assump- 
tion. 

We note that if one has the usual Kossakowski form 
of Markovian evolution [2], we can "diagonalise" the 
sets of operators to arrive at a Lindblad form [T3] where 
the decoherence operators are orthogonal and traceless. 
However, this standard form is not convenient for inver- 
sion, nor does it give much physical insight into the pos- 
sible processes leading to dephasing. The canonical form 
Eq.[l]decomposes the dephasing into operators represent- 
ing correlated level perturbations of orders 1 to N-l. 

V 4- CanonicDephasing (W) 

Calculate Canonical Dephasing Operators 

In: W Matrix (TV x K), feth column equals diagonal ele- 
ments of Lindblad operator Vk 

Out: V Lower triagonal matrix, columns equals diagonal el- 
ements of canonical Vk 



1 

2: 

3 
4 

5: 
6: 
7: 

8 
9 

10 
11 
12 
13 
14 
15 
16 

17: 
IS 

19 
20 
21 



R Number of rows of W 

C Number of columns of W 

V 4- W-ONES{R, 1) * W[l, :] 

k 1 // Running column index 

for r 4- 2,...,N 

i\ 4— Index 1st nonzero entry of V [r, k : 

il 4— ii + k — 1 // shift index 

while more than one element of V[r, k 



C] 
C] 



C] 



i-2 Index 2nd nonzero entry of V[r, k 
ii i'2 + k — 1 // shift index 
n 4- \V[r,h]\, r 2 4- \V[r,i 2 ]\ 

4>l 4- PHASE(V[r,ii]), (j>2 4- PHASE(V[r, i 2 ]) 



^V[:,i 2 ])/n c 
-^V[:,i 2 ])/n c 



n c 4- \Jr\ + r\ 

V[:,ii] 4- (rie +i4,2 V[:,h] + r 2 e^ 
V[:,i 2 ] 4- (rae-'^VMiJ-rie- 
if V[r,k : C] has non-zero entries 
fco 4— Index of 1st non-zero entry 
fco ko + k — 1 

V 4— Swap columns k and fco of V 
k 4- k + 1 

V <— Remove columns of V, apply phase corrections 



Algorithm 1: Canonical Dephasing Operators 
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